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The results of the searches for CP non-conservation in the decays of r leptons are presented. No evidence of violation of 
CP symmetry is observed neither in CLEO nor in BELLE data. Interpretation of these results is done within the framework 
of a model with a scalar boson exchange. Limits on the imaginary part of the coupling constant A, parameterizing the relative 
contribution of diagrams that would lead to CP violation, are -0.046 < 3(A) < 0.022 at 90% C.L. and 13(A) | < 1.7 at 90% C.L. 
in r — » irw°i> T and r — > K-K a v T modes, respectively. 



1 Introduction 

The violation of the combined symmetry of charge conju- 
gation and parity (CP) has been of long-standing interest 
as a possible source of the matter-antimatter asymmetry [ 
|l| in the Universe. Efforts to search for the CP-violating 
effects have concentrated so far on the hadronic sector. 
CP violation in strange meson decay has been the subject 
of intensive investigation since its first observation [ || in 
1964. Studies of hadronic [ [| as well as semileptonic [ 
[|, |j| |(| kaon decays provide precision measurements of 
the CP violation parameters. Searches for corresponding 
asymmetries in B meson decays are the focus of several 
large ongoing experiments [ f^, p|. Recent indications of 
possible neutrino oscillations [M make it important to 
re-examine the question of CP non-conservation in the 
leptonic sector. Such violation is forbidden in the Stan- 
dard Model but appears as a consequence of its various 
extensions [ [l0| . Among best-known theoretically are the 
multi-Higgs-doublet models (MHDM) [ (jj . Mod- 

els predicting lepton flavor violation often also predict CP 
violation in lepton decays J_OJl^]. Precision studies of 
muon decay parameters [ |16|, |17| show no indication for 
CP violation in such decay. 

In this talk I discuss three searches for CP violation 
in tau lepton decays. The searches for the CP-violating 
asymmetry in single tau decays were done by the CLEO [ 
[iJI and the BELLEf [l9) collaborations. In the third 
analysis by CLEO [ g(| correlated tau decays are used to 
search for a non-zero value of CP-odd optimal observable 
indicating CP violation. 



*Talk given at the Frontiers in Contemporary Physics-II (FCP01) 
Workshop, Vanderbilt University, Nashville TN, 5-10 March 2001. 



3350101-001 




Figure 1: Interference between two amplitudes with CP-even and 
CP-odd relative phases 8 and <f>. 

2 CP violation 

CP violation is the observable difference between a pro- 
cess (a — > b) and its CP conjugate (a — ► b). The neces- 
sary condition to observe such a difference is to have a 
non-zero CP violating phase <p in the Lagrangian describ- 
ing the process. Since the absolute value of a CP-odd 
phase does not have a physical meaning, CP violation 
can be observed only when interfered with a CP-even 
phase <5. Therefore, two interfering diagrams is needed 
to observe CP violation (see Fig. |l|). The probability 
density for such a process is given by: 

\A\ 2 = (Ai + A 2 e l<l 'e i5 )(Ai + A 2 e~ i<t 'e~ l5 ) 

= A\+ A\ + 2AiA 2 cos 4> cos 8 - 2A X A 2 sin sin (5. (1) 

The last, underlined, term is CP-odd since the phase 
<fi changes sign under CP conjugation. If this term is 
not equal to zero, physical observables between a process 
(a — > 6) and its CP conjugate (a — > b) may be different 
and CP is violated. 

A\ and A 2 denote the amplitudes and for physical 
processes they must be different from zero. Thus CP- 
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Figure 2: Amplitude for (a) standard W exchange, (b) scalar ex- 
change for t — > ■Kif'v-r decay. 



odd term is not equal to zero if the factors e 1 ^ and e tS 
differ from zero and are complex. 



3 CP violation in r decays 

A possible scenario for CP violation in r lepton decays is 
described [ |ll| by the interference of the Standard Model 
decay amplitude mediated by the W boson (amplitude 
Aw) with the amplitude mediated by the charged Higgs 
boson in the multi-Higgs-doublet model (amplitude Ah). 
These amplitudes play the roles of Ai and A2 of the pre- 
vious section (see Fig. |^). In this scenario, the charged 
Higgs couples to quarks and leptons with complex cou- 
pling constants and, thus, there can be a weak complex 
(CP-violating) phase (sin^> ^ 0). The overall Higgs cou- 
pling to the 7T7r° system is denoted by A: 

A = SR(A) + i9f(A) = |A|(cos^ + isin</>). (2) 

The usual choice for the CP-even phase 8 is a strong 
phase [ ^ which arises due to the QCD final state in- 
teractions between quarks. In the following, I consider 
only t decays into hadronic final states and a neutrino. 
There are several possible tau decay modes that can be 
used to search for CP violation. The r — > Kttv t mode 
is favored with the largest possible CP violation effects [ 
p^| . The corresponding branching fraction is, however, 
Cabbibo suppressed and the kaon identification efficiency 
in CLEO data is small, thus, leading to a small sample 
that can be used for experimental studies. This is not the 
case for r — > irir°v T decay due to the largest branching 
fraction of tau decaying into two pion state. However, 
scalar decaying into two pions implies isospin violation [ 
|22|. Therefore, the CP violating effects in this mode are 
suppressed, t In the following I discuss searches for CP 
violation in both tau decay modes. 



TThe suppression factor is proportional to the mass difference of u 
and d quarks. It is estimated to be of the order of 0.0035 [ E2]]. 
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Figure 3: Definition of angles in the hadronic rest frame. (Kir) 
defines the direction of the hadronic rest frame from the laboratory 
frame. 



4 Searches for CP violation in single tau decays 



There are two analyses [ [L8l 19 which use single tau 



decays into Kttv t and irir final state using the technique 
proposed by J.H. Kiihn and E. Mirkes [|23|, f22| . It can be 
shown that for the decays of unpolarized tau leptons the 
CP-odd part of the probability density is proportional 
to: 

P 0(W ~3(A)3(e l5 )cos/?cos?/>, (3) 

where 9(A) is an imaginary part of the overall Higgs 
coupling (see Eq. ||), 8 is a relative strong phase between 
W and Higgs exchange diagrams. /3 and ip are angles 
between a hadron, r lepton and direction of the hadronic 
rest frame in the laboratory frame (see Fig. |^). Therefore, 
the distributions of cos (3 cos ip for r + and t~ decays are 
given by: 



dN(r+) 
d cos (3 cos tp 



= const + ci9(A) cos /3 cos ip, (4) 



and 



dN(r~ 



d cos (3 cos ip 



= const + ci3(A*) cos /3 cos tp, (5) 



where const is contribution from the standard, CP-even 
part of the probability density and c\ is a proportionality 
constant. Asymmetry, defined as a difference between 
Eqs. |] and [| normalized by the sum of these distributions 
is: 



ci 
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f(A) cos/Jcos')/'- 
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Table 1: Asymmetry A cp for negative and positive values of 
cos /3 cos ip, CLEO Collaboration. 
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Figure 4: Asymmetry A cp as a function of cos/3cosi/> (CLEO). 



Non-zero slope in cos (3 cos ip distribution indicates CP 
violation. The value of 3(A) can be determined if em- 
pirical coefficient ci/ const is known.''' 



4-1 CLEO t — * K~kv t results 

The first search for CP violation in tau decays is per- 
formed in the r — > K~kv t mode [ Data collected 
from e + e~ collisions at a center of mass energy (y/s) of 
10.6 GeV with the CLEO II detector [ || at the Cornell 
Electron Storage Ring (CESR) is used in this analysis. 
The total integrated luminosity of the data sample is 4.8 
fb -1 , corresponding to the production of 4.4 x 10 6 tau 
pairs. A t~ — * K a s hTv T , K® — > 7r~7r + event sample is se- 
lected since three charged tracks in the final state are well 
measured. Here h~ is a charged pion or kaon. The stan- 
dard selection criteria is used to reconstruct r — > K®hv r 
by CLEO in previous analysis [ |25| . 

Observed asymmetry A cp (see Eq. ^) in the data 
sample for the positive and for the negative values of 
cos (3 cos ip after subtraction of the background is given 
in Table |l|. The distribution of A cp for the data and 
Standard Model Monte Carlo is illustrated in Fig. 

To relate the observed asymmetry A cp to the imag- 
inary part of the coupling constant A, the value of the 
proportionality coefficient c\ (see Eq. ^) is estimated us- 
ing modified TAUOLA package [ p6| . No CP violation is 
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tCan be estimated in Monte Carlo simulation. 



Figure 5: Asymmetry A C p as a function of cos/3cosi/j. The dashed 
lines are fits with a straight line or a constant(Belle). 



observed and the following upper limit on CP-violation 
is set: 

|9f(A)| < 1.7, at 90% C.L. (7) 
4.2 BELLE t -» vrvrV results 

There is an ongoing search [ |l9| by the BELLE Collab- 
oration for CP violation in r — > irir a v T mode using the 
same technique. The data sample used in this analysis 
has been collected in e _ e + collisions at a center of mass 
energy of 10.6 GeV with the Belle detector at the KEK 
asymmetric energy collider (KEKB) [|2?]. The total inte- 
grated luminosity accumulated during this period is 6.7 
fb -1 , corresponding to 6.2 million tau pairs. The final 
sample of r — > irir°v T contains 2.6 x 10 5 events. Using this 
sample the asymmetry from t + and t~ decays has been 
measured in two intervals of cos (3 cos ip similarly to the 
previous section. The values of the observed asymmetry 
are given in Table 0. The distribution of the asymmetry 
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Table 2: Asymmetry A cp for negative and positive values of 
cos /3 cos ip, Belle collaboration. 



cos (3 cos tp 


4 


cos (3 cos ip < 
cos (3 cos ^ > 


0.0056 ±0.0027 
-0.0004 ±0.0027 



is illustrated in Fig. ||. The value of the asymmetry is 
consistent with previous CLEO result and indicates no 
CP violation in r 
errors: 

| A* I < 0.010, at 90 % C.L.. 



7T7T v r decays within measurement 

(8) 



The limit on the imaginary part of the scalar coupling 
will be set in the near future. 



4-3 Summary of searches in single tau decays 

The first search for CP violation in r decays is done in 
t — * Knv T mode and provides for the first time a restric- 
tion on the imaginary part of a scalar coupling in tau 
decays: 

|9f(A)| < 1.7, at 90% C.L.. (9) 

Analysis using a single tau decaying into tttt°v t final state 
is in progress and current upper limit on the CP-violating 
asymmetry is: 



IA.pl < 0.010, at 90 % C.L. 



(10) 



This limit will be improved with better understanding of 
the systematics and an increase of the data sample. A 
conversion of the above upper limit to the imaginary part 
of the scalar coupling will be done in the near future. 

Searches for CP violation in a single r decay are sen- 
sitive only to a spin-independent CP-violating term in 
the decay rate. Other spin-dependent terms are propor- 
tional to tau polarization vector and, therefore, average 
to zero. However, the situation is different for decays 
of tau pairs produced in e + e~ annihilations, where the 
parent virtual photon introduces correlations of the r + 
and t~ spins. This correlation provides an additional 
constraint which permits to be sensitive to all CP-odd 
terms of the decay rate and maximizes the sensitivity to 
possible CP-violating effects. In the following, I discuss 
an analysis [ ||(| without these disadvantages. 

5 CP violation in correlated tau decays 

To maximize the sensitivity of our search an optimal vari- 
able is constructed, first proposed [ by D. Atwood and 
A. Soni, £, with the smallest associated statistical error. 




Figure 6: Correlated decay of tau pair into itk v t each. 

The variable is equal to the ratio of the CP-odd and CP- 
even parts of the total cross section assuming that the 
absolute value of the coupling A is unity: 

Podd 



p 



(ii) 



where P dd and P even are the CP-odd and CP-even parts 
of the total cross section. These terms are functions 
of the vector and scalar hadronic currents. A vector 
hadronic current is parameterized by the relative momen- 
tum between the charged and neutral pions multiplied by 
the vector form factor, /„, described by p Breit-Wigner 
shape: 



fv = 



(12) 



s — m 2 + imT(s) ' 
where s is a squared invariant mass of two pions, m n 
is a pion mass, and m and T(s) are the mass and the 
momentum-dependent width of the resonance, respec- 
tively. The latter is defined as: 

--^f' 2 ifs> (2m.) 2 



r( s ) 



o 



elsewhere. 



(13) 



Here, the contribution from the p' resonance [ ^9] is ne- 
glected. The scalar hadronic current is parameterized 
as a product of a dimensional quantity, M = 1 GeV/c 2 
providing overall normalization, and a scalar form factor 
f s . The choice of the scalar form factor is not unambigu- 
ous. Three possible cases are studied: one with f s = l, 
the second with f s described by the ao (980) Breit-Wigner 
shape, and the third with f s described by ao(1450) Breit- 
Wigner shape (see Eq. |l2|) with a width given by: 



r( s ) 



o 



(^,)V 2 if s >(2m^ 



elsewhere. 



(14) 



Since P dd is proportional to the imaginary part of A, 
then it can be expressed as: 



Podd(A) = %(A)P odd (l), 



(15) 
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> fs = 1 
■ O fs = BW[aO(9S0)] 
s = BW[aO(1450)] 



♦ tb = 1 

O f s = BW[a0(9B0)] 

• fs = BW[aO(1450)] 




Table 3: Average values of the optimal observable < f > for the 
Standard Model Monte Carlo and the proportionality coefficient c 
for the CP asymmetry fits for different scalar form factors. 



Form factor, f a 


< f >, l(T a 


c, 10~ a 


1 

BW(a (980)) 
BT^(a (1450)) 


0.7 ±0.6 
1.0 ± 1.1 
0.5 ±0.8 


66.8 ±4.3 
586.4 ± 19.4 
145.8 ±7.3 



Figure 7: Optimal observable £ for (a) Monte Carlo with no CP 
violation and (b) Monte Carlo with a maximal CP violation 3(A) = 
1. 



and 



< £ >= 3(A) f Po * d{l)2 dLips. 

J -Eleven 



(16) 



The integral in Eq. |l^ is always larger than or equal 
to zero, and equality occurs only if the odd part of the 
cross section vanishes. Therefore, the average value of £ 
is proportional to the imaginary part of the Higgs cou- 
pling constant and is positive if 53(A) > and negative if 
53(A) < 0. Monte Carlo simulation of the £ distributions 
for the three choices of the scalar form factors and no 
CP violation are shown in Fig. 0(a). The same distri- 
butions for the CP violating case 53(A) = 1 are shown 
in Fig. 0(b). The structure in these distributions is due 
to the resonant structure in the vector and scalar form 
factors. 

To relate the observed mean value of the optimal ob- 
servable < £ > to the imaginary part of the coupling con- 
stant A, the 53(A) dependence of < £ > must be known. 
To the first order the mean value is proportional to 53(A) 
with a proportionality coefficient c: 



< £ >= c x 3(A). 



(17) 



The value of this coefficient is estimated using Monte 
Carlo simulation as described in the Section 5.2. 



5.1 Data and Monte Carlo Samples 

The data used in this analysis were collected at the Cor- 
nell Electron Storage Ring (CESR) at or near the energy 
of the T(4S I ). The data correspond to a total integrated 
luminosity of 13.3 fb _1 and contain 12.2 million t + t~ 
pairs. Versions of the CLEO detector employed here are 
described in Refs. [ |3| ||o|]. From this data sample, events 
consistent with e + e~ — > t + t~ interactions are selected 
where each r decays into the tttt ^ final state. The 
event selection follows mostly the procedure developed 
originally for the study of the r — > pv T decays [ In 
the following section it is shown that the event selection 
criteria do not introduce artificial CP-violating effects. 



To estimate backgrounds large samples of Monte 
Carlo events are analyzed following the same procedures 
that are applied to the actual CLEO data. All non-tau 
backgrounds are found to be negligible. The main re- 
maining background is due to the r-pair events in which 
one of the r's decays into pv T while the other decays into 
7r27r°^ r and the photons from one of the ir°'s are not de- 
tected. The contamination from this background source 
is estimated to be 5.2%. The second largest background 
contribution of 2.1% is due to one of the r's decaying 
into the K*v r final state producing a neutral pion plus 
a charged kaon which is mistaken for a pion. All other 
tau decays provide much smaller contributions with the 
largest being less than 0.7%. The total background con- 
tamination from tau decays is estimated to be 9.9%. This 
background neither introduce a bias in the < £ > distri- 
bution nor modifies the value of the coefficient c, defined 
in Eq. [P?]. 

5.2 Calibration 

To relate the observed mean value of the optimal observ- 
able < £ > to the imaginary part of the coupling constant 
A, coefficient c must be determined. It is estimated us- 
ing several signal Monte Carlo samples generated § with 
different values of 53(A). For each sample we calculate 
the average value of the optimal observable £ and plot it 
as a function of 53(A). For each form of the scalar com- 
ponent, the calculated asymmetry distribution is fit to 
a straight line to obtain the calibration coefficients. To 
check that the selection criteria do not create an artificial 
asymmetry the mean value of the optimal observable is 
calculated for Standard Model Monte Carlo samples for 
each choice of the scalar form factor. These values along 
with calculated coefficients c are listed in Table |[ For all 
three form factors, the mean value of £ for the Standard 
Model Monte Carlo sample is consistent with zero within 
its statistical error. Therefore, event selection criteria do 
not introduce an artificial asymmetry in the £ distribu- 
tion. The coefficients from Table || are used to calculate 
3(A). 



s These events are generated with full GEANT- based detector sim- 
ulation [ re2| and pattern recognition software. 
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Figure 8: The distribution of the CP violation sensitive variable £ 
for the data (dots) compared to the Standard Model Monte Carlo 
prediction (solid line) for (a) / a =l, (b) f s = BW(a (980)) and (c) 
f s =£W(a (1450)). 



Table 4: Average value of the optimal observable § after subtract- 
ing the average value for Standard Model Monte Carlo, calculated 
value of 9(A) and 90% C.L. intervals on Of (A). 



/. 


< f >, io~ a 


3(A), 10-* 


3(A), 90% C.L. 


1 

a (980) 
a (1450) 


-0.8 ±1.4 
-0.6 ±2.4 
0.2 ± 1.7 


-1.2 ±2.1 
-0.1 ±0.4 
0.1 ± 1.2 


(-0.046; 0.022) 
(-0.008; 0.006) 
(-0.019; 0.021) 



5.4 Summary 

No sizable systematic errors are found\diich could alter 
the limits shown in Table ^. Within our experimental 
precision no significant asymmetry of the optimal vari- 
able is observed and, therefore, no CP violation in t de- 
cays. Due to the uncertainty in the choice of the scalar 
form factor the most conservative 90% confidence limits 
are used corresponding to f s = 1: 

-0.046 < 3(A) < 0.022, at 90% C.L.. 

These limits include the effects of possible systematic er- 
rors. 

5.5 Pseudo-helicity method 

In this section I discuss another method to search for 
CP violation in r decays. It is more intuitive but less 
sensitive than the optimal observable method described 
above. 

The helicity angle, 6^, is defined as the angle be- 
tween the direction of the charged pion in the irir rest 
frame and the direction of the irir system in the r rest 
frame. In Standard Model, the helicity angle is ex- 
pected to have a distribution corresponding to a vector 
exchange: 

— a + focos 2 ^. (18) 

a cosUktt 

For scalar-mediated decays, there is an additional term 
proportional to cos 9^ that corresponds to the S-P wave 
interference and linearly proportional to the scalar cou- 
pling constant A. In general, A is complex and the term 
linear in cos 9^ is proportional to the real and imaginary 
parts of the scalar coupling with coefficients c\ and C2, 
respectively: 



5.3 Observed 



values 



For each choice of the scalar form factor, a distribution 
in £ is obtained. These distributions are shown in Fig. ||, 
with those from the Standard Model Monte Carlo simula- 
tions overlaid. From these distributions the mean values 
< £ > are computed after subtracting the average value 
for Standard Model Monte Carlo, which are reported in 
the first column of Table ||. In each case, the appropri- 
ate empirically-determined coefficient given in Table ^ is 
used to derive a value for the imaginary part of the Higgs 
coupling A, as described in the preceding section. These 
values along with the 90% confidence limits on 3(A) are 
reported in the second and third columns of Table 0. 



; — ^ — ~ a+ci!ft(A) cos6> 7r7r +c23(A) cos 0^+6 cos 2 0^. 
d cosO-x-K 

. (19) 

The observation of the terms proportional to cosine of 
the helicity angle would indicate the scalar exchange in 
the tau decays. 

In order to calculate the helicity angle tau rest frame 
must be known. Due to the unobserved neutrino, the tau 
rest frame can only be reconstructed with a two-fold am- 
biguity. Such ambiguity can be avoided by using the 
pseudo-helicity angle, 9*. This pseudo-helicity angle is 
obtained by replacing the tau rest frame with the labo- 
ratory rest frame where it is defined as an angle between 
the direction of in the tttt rest frame and the direc- 
tion of the tttt system in the lab frame. The difference 
between the pseudo-helicity distributions for the r + and 



^The detailed description of the most significant effects is given 
elsewhere [ BOf . 
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Figure 9: (a) Pseudo-helicity distribution for r~ and r + in data, 
(b) difference between pseudo-helicity distributions for r~ and r+ 
in data. Solid line is a linear fit and dashed lines show extrapolation 
to the region excluded from the fit. 



dependence as a function of 3(A) is fit to a straight line 
to obtain c 2 : 



c 2 = (107.5 ± 12.6) x 10" 



(21) 



This coefficient is used to obtain the value of the imagi- 
nary part of the scalar coupling 3(A) after the correction 
of the slope by Monte Carlo: 



3(A) = 0.028 ±0.037, 



(22) 



and 



0.033 < 3(A) < 0.089 at 90% C.L. (23) 



As expected, the limit on the 3(A) is less strict than the 
one obtained using the optimal observable. 



t~ decays is expected to have the same form as given by 
Eq. [l^ but with a different numerical coefficients. 

The term including 3(A) changes sign for tau lep- 
tons of opposite charges. Therefore, the difference of 
the pseudo-helicity distributions for positive and nega- 
tive tau leptons has the term linear in cosf?* proportional 
to the imaginary part of the scalar coupling A only: 



dN(r~) dN{r + ) 



2c 2 3(A) cos6>*. 



(20) 



The presence of this term indicates CP violation. 

In this study the same data sample is used as for 
the optimal observable analysis with the same selection 
criteria. The pseudo-helicity distribution for r~ and r + 
is given in Fig. ||(a). 

The structure in Fig. ^(a) is due to the variation in 
the efficiency as a function of charged pion momentum 
and 7T° energy. To obtain the product of the imaginary 
part of the scalar coupling A and a linearity coefficient 
C2 (see Eq. |o|) , the difference of the two pseudo-helicity 
distributions for negative and positive tau leptons is fit to 
a first order polynomial. To minimize systematic effects 
due to soft pion reconstruction the fit is performed in the 
region of —0.7 < cos 9* < 0.8, which corresponds to pions 
with momentum higher than 0.3 GeV/c. The obtained 
value of the slope for the data distribution is: 

c 2 3(A) = (4.2 ±3.6) x 10~ 4 

and illustrated in Fig. ||(b). It is consistent with zero 
within statistical error. 

To determine the proport iona lity coefficient c 2 the 
procedure described in Section 5.2 is followed. The slope 
of the difference of pseudo-helicity distributions is calcu- 
lated after applying selection criteria. Five samples of 
200,000 signal Monte Carlo events generated with differ- 
ent values of 3(A) are used. Then, the calculated slope 



6 Summary 

I have discussed three searches for CP violation in the 
decays of tau leptons. No evidence for CP violation is 
observed within the experimental precision. First two 
analyses search for non-zero CP-violating asymmetry in 
the single tau decays to the K-kv t and ttk v t final states. 
The following upper limit on the imaginary part of the 
MHDM scalar coupling constant A is set in r —* Kirv T 
decays: 

|3(A)| < 1.7 at 90% C.L, (24) 



The value of the asymmetry for r — > irir°i/ T is: 
\A cp \ < 0.010, at 90 % C.L. 



(25) 



The conversion of this asymmetry to the value of 3(A) 
will be done in the near future. 

The third analysis used optimal variable to search 
for CP violation in correlated tau decays decaying into 
irK°v T final state each. The following confidence limits 
on 3(A) are set: 



- 0.046 < 3(A) < 0.022 at 90% C.L. 



(26) 



All these results agree with each other and the limits 
on imaginary part of the coupling constant A restrict 
the size of the contribution of multi-Higgs-doublet model 
diagrams to the r lepton decay. 
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